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Abstract 



Supersymmetric standard model has a parity violation in QCD through 
chiral quark-squark-gluino interactions with non-degenerate masses be- 
tween left-handed and right-handed squarks. Since experiments have not 
show any parity violations in QCD yet, a bound for the mass degeneracy 
between left-handed and right-handed squarks should exist. In this pa- 
per we try to obtain this bound for each squark. At first, we investigate 
a non-degeneracy bound between m^ L and m 5H from experimental data 
of charmonium decay. Next, we estimate the non-degeneracy bounds 
for u and d from nucleon-meson scattering data, and comment on other 
squarks. 



1 Introduction 



Supersymmetric (SUSY) standard model (SM) is one of the most promising candidates beyond 
the SM due to a gauge coupling unification, an existence of dark matter, and so on. It is worth 
noting that the SUSY SM has chiral gauge interactions in fermion-sfermion-gaugino vertexes, 
for example, left-handed fermion only couples to left-handed sfermion with gauginos. We know 
this left-right index in sfermion is just a label and has nothing to do with spin, but anyhow, 
fermion-sfermion-gaugino interactions are exactly chiral. Moreover, a mass of left-handed 
sfermion is not the same as that of right-handed sfermion in general, and this non-degeneracy 
is also induced from radiative corrections picking up weak interaction. Therefore, due to 
this non- degeneracy, SUSY gauge interactions cause parity violation even in QCD! This is a 
remarkable feature of the SUSY SM, and this effect is negligible in a lot of other candidates 
of beyond the SMpQ. We take a setup of i?-parity conservation, where spartcles propagate 
only inside loop diagrams due to their heavy masses. Since any parity violations have not 
been discovered in QCD, non-degeneracy bounds should be obtained between the left- and 
right-handed squark masses. Does an experimental fact of well conservation of parity in QCD 
suggest a good degeneracy between left- and right-handed squarks? We must investigate this 
degeneracy, and try to obtain the non-degeneracy bounds from current experimental data. As 
for the degeneracy of II and Ir, there are some researches in a process of tt pair production 
in collider experiments [U [2]. Where amplitudes of tt pair production depend on heli cities of 
them, and non-degeneracy between ti and tR causes the asymmetry measurement of the cross 
section. A case of 0(100) GeV (0(1) TeV) masses of sparticles was investigated in Ref.[2] 

(Hrf.H!). 

Here we should comment on sparticle masses which are consistent with collider experiments. 
One case is that all sparticles are heavy of 0(1) TeV as well as gluino mass > 700 GeV[3], 
where all sparticles are too heavy to be detected at detectors in the current experiments. 
Another case is that light sparticles exist which are degenerate within 30 GeV other than the 
heavy sparticles. It is because there are experimental cut for pts of multi-jets with missing 
transverse momentum in SUSY search at LHC (Tevatron), where an event selection for jets is 
p T > 40 GeV[3] (j>t > 30 GeV[l]), and pt of jets are roughly estimated as the mass difference 
of gluino and squarks. Thus, the degeneracy of sparticles within 30 GeV is consistent with 
collider experiments. We will take both two cases in the following numerical analyses. 

In this paper, we try to obtain the bound for left-right degeneracy of squark masses other 
than stop. At first, we investigate a non- degeneracy bound between mg £ and m £fl from exper- 
imental data of charmonium decay. For this analyses, we use non-relativistic QCD (NRQCD) 
[5], since the charmonium is enough heavy to use it. We also estimate the non- degeneracy 
bounds for u and d from nucleon-meson scattering data, and finally comment on bounds for 
other squarks. 
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2 Dimension six operators from SUSY 

We calculate dimension six operators by intgrating out sparticles in the framework i?-parity 
conservation. By integrating out SUSY particles, we can obtain higher order gauge invariant 
operators in terms of the SM fields. We calculate them up to 0(a 2 s ) and neglect 0(a s a y ), 
where a s (a y ) is gl/^K (y 2 /in) with a QCD (Yukawa) coupling, g s (y). It is because Yukawa 
couplings, ys, which we deal with are all smal£], and up to this order, sfermion left-right 
mixings are negligible. In the NRQCD, there are three types of dimension six operators, 
Of Fl and OqgG, which can contribute the parity violation in QCD. Here O^p and Ofp 
are color singlet and octet 4-Fermi operator, respectively. Other dimension six operators such 
as q-q-G-G and q-q-G-G-G vertexes are next leading order in the NRQCD, so we neglect them 
in the following discussions. 

For O^p and Of F , they are given by [6] 

„ 4 flavor 

O$0«0 = ^ E ^ (^ p lQ) tiltPuf) + Crr (qi^Prq) ($VW) 

+C LR {q-fP L q) (q'^Pnq) + C RL (q^Pnq) (qi,P L q)] , (2.1) 



flavor 



= ^jjh E ^ Dll (^V^l<?) (^V^) + Drr (qT a YP R q) (q'T a ltl P R q') 



+D LR (qT a ^P L q) tfT a ^P R q') + D RL (qT a ^P R q) tfT a lfi P L q')} , (2.2) 
where coefficients are 

2 

Cll = g[fi(mg L ,rrig0 + / 2 (m & , m 9j )\, (2.3) 
2 

C RR = g\fi{mq R ,mg> R ) + / 2 (m fe , m^)j, (2.4) 
2 

C LR = -Q\fi(m 9B ,m ¥x ) - f 2 {mq L ,mq /R )}, (2.5) 
2 

C RL = —g[fi{m 9L ,m? R ) + / 2 (m & , m^)], (2.6) 

1 7 

£>ll = -g/i(m & ,m5j.) - gMmfe,m^), (2.7) 

1 7 
£>i?i? = - g /1 "V H ) - g/2(m^,m^), (2.8) 

7 1 

^Lfl = --/i(m g - L ,m^) - -/ 2 (m^,m^), (2.9) 

7 1 

D RL = --/ 2 (m 9 - H ,m^) - ^fi{mg R ,mg' L ), (2.10) 



We do not analyze t in this paper. 
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f 1 f 1 y Z 2 

fi(mq,m?) = dy dz— Y — — , (2.11) 

Jo Jo {m z ~-mf)yz+(m z ~-m z ~,)z + rriq, 



"g q/u • \ q q' 

/2 (m 5 , m? ) = J g dy^ "z [{ml _ mf)yz ^_ ml)z + mf]2 . (2.12) 
On the other hand, O qq c is given by 

(2.13) 

= (e^m^)^! + e 2 {m q - L )$2}K + {ei(m^)^ 2 + 62(771^)^1}^ 

+ {e 3 (m g -J(fc? + fcf) - e^mqjh ■ k 2 }^ - e 5 (m 9 ~Jze afttl/ 7 5 7^ia£;2/3, (2.14) 
E» = E»(m qR ), (2.15) 

ei(m^) 

107m? - 495m|m? + 477m?m| - 89m? - 6(m? + 3m~m~ - 54m?m| + 18m?) log(m?/m~) 
" 18 (m? - m?) 4 ' 

(2.16) 

e 2 (777g) 

—203m? + 351m|m? — 189m?m| + 41m? + 6(m? + 51m~m~ — 54m?m| + 18m?) log(m?/m~) 
" 18(m? -m?) 4 ' 

(2.17) 

e3^) = e 2 (mg), (2.18) 
e 4 (m,j) 

155m? + 423m~m~ — 333m~m~ + 65m? + 6(m? + 27m~m? — 54m~m~ + 18m?) log(m?/m^ 



9(m? - m?) 4 
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(2.19) 

e = 9(m| - m| - 2m|m? log(mf /mf )) ^ q 

(mj — m~) 3 



We will use these QCD dimension six operators in order to obtain the non- degeneracy bounds 
of left-right squark masses. 



3 Quarkonium in NRQCD 

Let us formulate a quarkonium of qq bound state in the NRQCD framework by introducing 
a bilocal field. It is applicable for heavy quarks, and a related work has been shown in, for 
example, Refs. El U Ei] • 
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At first, we integrate out gluon field, then, QCD action becomes 

S= [ iq(i$-m)q] + (-i) [ f jI(x)V^(x - y)j v (y), 



(3.1; 



x Jy 



where T>^ u (x — y) is gluon propagator. When q is a heavy quark, it is expanded by its mass as 



q(x) 



ipe- imt + %^xz imt 



(3.2) 



(f and x denote particle and anti-particle components, respectively, and this expansion is 
so-called Foldy-Wouthuysen-Tani transformation [TTj IT2"] . In non-relativistic limit, the gluon 
propagator in the second term induces a (gluon) potential as 



Ty{x - y) 



(2vr) 4 ^-b1 



An\x — y\ 



(3.3) 



Taking a color singlet part in the second term of Eq. (13. ip (color octet part is the next leading 
order [5]), we can obtain a NRQCD action, 



Snrqcd 
1 



2m 



^ f ido + 7T- J <P + X 1 f id - — 



2m 



X 



+ 



2N. 



c 




5(x° - y°) [p\x) X {y)V{r) X \yMx) + V \x)a\{y)V{r)x\y)a^{x)\ , (3.4) 



x Jy 



where r = \x — y\ and V(r) = C7p^/(47rr), and color factor comes from T£Tj£ = -^f^biiSkj — 
-j^T-jTfoj through Fierz transformation. Hereafter, we note Nc as a color number, which is, 
of cause, Nc = 3. Next, by inserting the following identities, 



1= I l[Vs»Vr j expz I 1 4>l(x,y)(s^x,y)-^(x)^x(y)), 
l[Vs^Vr exp* I I ^(x,y)(s^(x,y)-xHx)a^(y))i 

.. ,, J x J y 



(3.5) 
(3.6) 



into Eq.( l3.4p . the QCD action becomes 
Snrqcd = 




x Jy 



[(p ] (x)K w ^(y) + x\x)K xx x(y) - 4>^(x, y)^ ] \x)a M x(y) - X^i^a^y)^!, y) 

(3.7) 



+ ^T 6 ( X ° ~ y°) s ^( x ' y)V(r) SlI {x, y) + 0^(x, y)s IM (x, y) + s j {x, y)^{x, y) 



where the kinetic terms denote 



V 2 \ 



(p\x)5 4 (x - y) \ id + —J ip(y) = ^ ] (x)K w ^(y), 
X*(x)6*(x - y) (id - x{v) = xH x ) K xxX(y)- 



(3.8) 
(3.9) 
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An effective action of the bilocal auxiliary field ^(x, y) will be obtained by integrating out 
s M , ip, and x- A potential term is induced by integrating out s M as 

P\x,y) [2N c 5(x° -y°)V-\r)g, u ] cj?{x,y). (3.10) 




x Jy 



On the other hand, tp- and ^-integrations will derive a kinetic term of M as shown below. By 
integrating out ip and \ i n Eq. fl3.7D . we can obtain the term 

iirl ° g V K X9 K xx ) ~ Zirl ° g { K xx ) + ^ n Ll { K£K W ) ' 

(3.11) 



where K vx and K vx are denoted as K vx = —(p^(x,y)a^ and i^ x</3 = —a^i^x.y). The Trlog 
is expanded in Eq. fl3.lip . where n — 1 is vanished by a trace, and the leading term is coming 
from n = 2. After taking traces of spinor, color, and coordinate indices, the leading term in 
Eq. fl3.11l) becomes 

Zl Tr ( K ~,l K ^ K xl K x, 
9 \ K~ X K K~ X K 



-iN, 



c 



/ / / / Tls p™ K w( x iy} K vx(yi z ) K xx( z ' w ") K xA w i x )' ( 3 - 12 ) 

JxJyJzJw 



where propagators are given by 



K -l( x y)= f Al l -ip{x-y) 5 «p 

v) = -f f2 d l o L- e^P*. (3.13) 



We use a center of mass coordinate and relative coordinate (0,r) M as x M = X M + |(0, r)'* 
and |/ M = A M — |(0, r)' 4 . The relative coordinate does not have tzme-component, since (p^(x, y) 
is a coincident bilocal field for x and y. Then, ^(x, y) is represented by 

= { Ul{l)e- lkX e- a ^ = [ [ [ttfle-nxe** (3.14) 

Jk J i Jk Ji° JT 

with their momentums as = (^- + 1°, f + and = (^ — i°, | — Z). In this frame, Eq.f l3.12p 
is written as 

-2iN c [ !J f— -1 r€\r)M^)e-^ +g) (3.15) 



'kJiJfJs k° — -r 

4m m 
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by integrating 1°. Then, the effective action of ^ is given by 



S, 



eff 




>X 



V(r) 



K x (r) 



^x{ 



(3.16) 



where Kx(r) = id x — ^ L — — . We omit overall factor 2Nc by use of normalization of the 
field. Note that a Green function (4 > xi^) ( t > Y (*)) i s given by 

^^t(s)) = [V^ 1 - K-\l (X, f; Y, 5) 

= V{r) 9llu 5\X - Y)S 3 {r- s) + [V(K - V)~ l V] (X, r; F, s). 



(3.17) 



In asymptotic states, X ^ Y, the first term vanishes. The second term is what we want, and 
V is rotated out by field redefinition, then Eq. fl3.17p becomes 



«W)=[(K-V)- 1 ] ia ,(X,?;Y,g). 
This means that the effective action in Eq.f l3.16l) can be rewritten as 



S, 



cff 




^(f)[K x (r)-V(r)}^ x (f). 



(3.18) 



(3.19) 



x Jf 



We should notice that this form is correct when asymptotic states exist [13] and <^x(0 is an 
on-shell state. This is the situation we consider hereafter. 

Now we estimate spectra of bound states. Since the potential in Eq. (13. 19ft only depends 
on relative coordinate, we can divide a motion of the bound state between that of relative 
coordinate and that of center of mass coordinate. By defining H(r) = — V^/m — V(r), the 
effective action in Eq. fl3.19p is written as 



S, 



cff 




X Jf 



id x -^JL + H(r) 



4m 



4>vx(r). 



i M can be expanded by a complete set of relative coordinate of ip n (r) as 

d 3 P 



J x^ 



E^(x)w l (o = E / 



(2k) 



r <(P)^ n (r)e- 



PX 



(3.20) 



(3.21) 



where a%(X) is a plane wave, and ip n (r) is a possible bound state which this system can take. 
An eigenstate of H(r) which satisfies 



H(r)%l) n (r) = E n tp n (r) 



(3.22) 



is a quarkonium, and E n denotes a binding energy of it. Orthogonality and completness 
suggest 



d 3 ril)i(f)ij m (r) = 5 nm , ^ ^n(r)^(s) = 5(f- s). 



(3.23) 
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A hadron wave function is factorized by a£(X), which only depends on center of mass coordi- 
nate. Here /i represents spin singlet (triplet) state of meson when fi = (fi = i). Note that a 
hadron labeled by n is created by a$(X) as a$(X)\Q) = \n). 

Here let us apply this formalism to a charmonium, for example. We denote n = rj c , h c , J/ip, 
Xci ' ' ' , then a spin singlet state ^{r) and a spin triplet state (fixi^) are represented by 

<&(f) = <(X)W(r) + al(X)i; hc (r) + • ■ ■ , (3.24) 
<PW) = a^(A> J/v ,(r) + a^jX^if) + ■■■, (3.25) 

respectively. We now obtain the effective action of charmonium in the SM QCD, where parity 
is conserved. 



4 Direct parity violation 

In the SUSY SM, parity can be violated in a quarkonium through the non- degeneracy of 
left-right squark masses. As we have shown, there are three parity- violating operators, O^p, 
Ofp, and O qq G- At a direct decay vertex of quarkonium, O^p gives the leading order of parity 
violation, and we call this process "direct parity violation" . The explicit form of the direct 
parity violating operator is given by 

O p 4 ; = (A uc + B cu )5\x - y)[u(x) lfl u(y)][c(x)Yl 5 c(y)], (4.1) 

where A uc and B cu are 

A 1 f „(u,g) . ^((«,c) , ^(u,c) ^-y(m,c)\ (a o\ 

= 1927T 2 4 RR LR R.L )' 

R jjgg 1 i n {c,u) n {c,ii) n {c,u) n {c,u)^ (A Q\ 

= 192?r 2 4^~ LL RR ~ L lr + c rl J> l 4 -^ 

respectively. We estimate -u-quark contribution at first, and later include ci-quark contribution. 
Note that squark flavor is labeled by C\f = L,R), and has squark mass dependence 

through fi(mq,mq') and /^("^g, Ttg')- For example, C^/ is denoted as 

2 

C ll = g [h {mu L ,mc L ) + f 2 {mu L , mgj] , (4.4) 
and other C-factors are similarly obtained by using Eqs. (l2.3pMI2.6p . 

(8) 

As for OggG and 0\p, they do not induce the leading order contributions, because they 
must emit a gluon in the decay vertex. We can neglect gluon exchange between in-going and 
out-going states at the decay instant in the NRQCD, since non-relativistic bound states are 
hadronized by space-like gluon exchanges. Therefore, we can neglect the contributions from 

(8) 

OqqG and 0\p, and factorize this decay process by a vacuum insertion as in Fig. 1. 
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time 




Figure 1: Factorization by a vacuum insertion in direct parity violation process. A box stands for O^p 
Eq.gU). 



We focus on a charmonium, r) c , which is h under J PC , and has mass of 2980 MeV. Notice 
that 0\p is a contact interaction, where the decay constant is a value of wave function at an 
origin due to 5-function and a decay through the contact interaction is only possible with the 
S-state (angular momentum L — 0). Thus, reminding tt is 0~ + , t] c can not decay to tttt until 
it pick up parity violation, since ir(p), ir(—p) systenill of S-state is ++ . Note that there exits 
weak interaction, however, it also breaks C. Anyhow, as in Fig. 1, the direct parity violation 
through the SUSY effects, i.e., a two-body decay process, r\ c ixix, should be factorized as 
{7m\Olp\r] c ) ~ (7r7r|97 M Q'|0)(0|g7 /J 7 5 g|?7 c ). Here (7T7r|g7 M g|0) is a pion form factor, and we can 
estimate (0|g7 At 7 5 g|?7 c ) by use of NRQCD. Actually, by regarding g"7°7 5 g ~ — |xV + h- c - i n a 
non-relativistic picture, the S-matrix element of r] c — > Tin is given by 

(nn\Oi;\ Vc ) ~ - 1 -{A UC + B cu )5\x - y)(7r7r\u\x)u(y)\0)(0\ X \x)cp(y)\r} c ). (4.5) 

Here (irir\vt(x)u(y)\0) ~ F s (k) is a scalar form factor of pion, which has non-trivial energy 
dependence. 

In general, when a bound state |n)(= a^(P)\0)) decays through a bilocal operator C^ A "(r) = 
0^(r*)0^(r*) • • • , its matrix element is given by 

(0\T[Of(mn) = d 3 sF-(Y; s) (idyo - S _ H( s )\ (0\T[O$"'(?)<# {a)]\0), 

where Fp(Y; s) = ip n (s)e~ tP ' Y , and it satisfies J d 3 Xd 3 rcj) x (r)Fp(X;r) = aJj(-P) from orthog- 
onality and completeness. Thus, the transition amplitude in Eq. fl4.6p is given by ip n {r)e~ lPX 
with O u x x -(r) = <p x (r). 

t It has P = (-l) L and C = (-l) s+i 
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Let us go back to a charmonium, and take q as c-quark in Eq. fl3.2p . Since a heavy quark is 
non-relativistically expanded as Eq.f l3.2p . the 4- Fermi operator can be also expanded similarly. 
In the leading order, components of x f an d <PX i n the bilocal field, are only creating and 
annihilating operators of charmonium. Thus, 4>xi^) corresponds to X^{ x )^{y)i an d we name 
a label n = r] c for the charmonium, which suggests 

<0|&(f)|»fc> =^e~ ip - x . (4.6) 

Remind that O^p is a contact interaction, and we can use m nc for an energy of the pion form 
factor due to a momentum conservation. Then, we obtain 

(ixixPi;^) ~ -^(A uc + B cu )F s (m Vc )iJ Vc (0)- (4-7) 

There is a d-quark contribution as well as u-quark ones, so that the effective 4- Fermi operator 
0\p becomes a linear combination of u and d. Therefore, T(r] c — >■ nn) is estimated as 

TY„ -nr^ 14 A +R R ,2 I^KJ| 2 |^ c (0)| 2 ( . 

J- \J]c ~+ ttttJ ~ A uc + A dc + # cu + i^cd — — . (4.8) 

lomi 

Note again that a decay through the contact interaction is only possible with the S-state, and 
in states with non-zero angular momentum. The wave-function at origin vanishes, and decay 
constant also vanishes. This is a characteristic feature in the direct parity violating process in 
the SUSY SM. 



5 Indirect parity violation 

The QCD dimension six operators from the SUSY SM can have the parity violating effects, 
and actually, they can also contribute organization of quarkoniums themselves. We call this 
effect "indirect parity violation" , and we investigate it in this section. For this indirect parity 
violation, all O qq c, O^p, and £)fp contribute as in Fig. 2. 

The indirect parity violation induces a mixing between an even-parity state and an odd- 
parity state as well as a S-state and a P-state in a quarkonium. As the parity violating term 
is written by SV^ u (r) in the potential, the effective action in Eq. fl3.19p includes indirect parity 
violation by rewriting V{r)g^ u — > V(r)g^ u + 5V tJ,v {r). Here 5V^ u (r) is a matrix in a basis of 
S- and P-states, which has off-diagonal elements of hadron state labeled by n (and /i). Now 
let us calculate the mixing between asymptotic states in the SUSY SM by using the basis of 
the SM QCD. Since the potential only depends on relative coordinate, the wave function can 
be expanded by ^ n {r) in the SUSY SM as 

r x (r) = J2 A n(^n(r), (5.1) 
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Figure 2: Diagrams which contribute indirect parity violation through dimension six operators (which are 
shown a box at a vertex). (Left): a contribution from O qq c, (Right): a contribution from O^p or ofp. 



where ^ n (r) satisfies eignvalue equations, 

[#QCD (f) + sv(7)}y n (f») = E^ n (f), (5.2) 

for E^ 1 ^ E n . Note that n is the label of the hadron, which contains an information of spin 
(/x = 0: singlet, \i = i: triplet). This ^ n (f) must be \I/ n (r) — > ip n (f) as 5V — > 0, so that it is 
given by 

^n(r) = Vn(r) + Yl 1F^V^*(0, (5-3) 
up to the first order of perturbation. Note that \I/ n (f) must satisfy 

J d 3 r^ n (f)^ m (r) = S mn (5.4) 
for the zeroth order of perturbation. V nk is defined by 

V nk = ! d 3 s4(s)[5V(s)}Ms). (5.5) 



The classical complete set \l/ n (r) should be written by the QCD complete set t/) n (r), while a 
hadron creation operator is given by Af£(X), so that A£(X) corresponds to a£(X). ^ n (0 and 
if] n {r) are different complete bases as 

= E<w^w = ( 5 - 6 ) 

Thus, by use of orthogonalization of *ff n (r), we obtain 

^w=< f w+E (v ri ( y • ^ 
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Let us consider a charmonium system. Equation (j5.7p means an observed t] c is almost 
represented by a mixing state of r\ c and Xco as 



l»7c) 



obs. 



\Vc) + 



VcXoO 



E - E 



|Xco)- 



(5.* 



Xco is ++ with mass of 3415 MeV, and a decay of Xco to 7r7r is possible (see, Eq. fl6.2p ) when 
7T-7T system has angular momentum, L = 1. We estimate parity violating potential induced 
from the SUSY SM. As for O^p in Fig. 2, its coefficient only depends on and rric R , since 
the bound state is charmonium. The parity violating terms in Ofp are given by 



192tt 2 2 



(5.9) 



where we use spin relation, 6 a p6 y \ 



i< A a^, and aV 6 



5 



a6 



ie abc a c 



factor is rewritten as \5ij5u 



A color 



2N C 



bu&kj + Tfi-TA for an exchange of spin. We must be careful 



for exchanges of spin and coordinate, where only spin-singlet changes its sign (Table{TJ). After 





exchange of spin (ip -h- x) 


exchange of coordinate (x -H- y) 


spin singlet 4>°(x,y) 
spin triplet (j> l (x,y) 


asym. 
sym. 


sym. 
asym. 



Table 1: exchanges of spin or coordinate 



careful calculations, O^p is given by 



-»• 



12^ 1 



192vr 2 2 



t 



(5,6) 
LL 



(c,c)> 







4^;V(r) 



4m c N c 

4V(r)$ 
4ze ij ' fc VrV(r 



(5.10) 



where V 



y) and 0* 



As for we can use the calculation result of 



0\p, since spin structure is the same. The different point is just color factor, and by using 
TATA = ^5 a 5 kj -^T^T^ we show color octet part is C F (= (N%-1)/(2N C )) times larger 



■ ij kl 



than dQ. Then, non-relativistic potential from 0± F and Ofp with parity violation is totally 
given by 



12g 



192tt 2 8m c JV c 



-C 



(c,c) 
LL 



c 



(c,c)> 
RR > 



LL 



'RR 





4fr'V(r) 



4ie^ fe *§~*V(r) 



(5.11) 
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For a non-relativistic potential from O qq c in Eq. fl2.14p . we estimate leading part. Since 
O qq G is not the contact interaction as O^p, its parity violation effects should be added to the 
gluon potential. The bilocal operator after integrating out gluon is given by 

C ~ W)^&l*Plm(*)] V(r) [q{y)T A l°q{y)\ , (5.12) 

where E° L R has eight terms in total, which are categorized as 

(i) ± ei {m,)-^— [fol^fl^fe) + (d»d°q x T A lfl q x )} V(r) [q(y)T A l°q(y)} , (5.13) 

(ii) ± e 2 ( m? )^|L_ [(d»q x T A ^d°q x ) + (d%T A lfl d»q x )} V(r) [q(y)T A l°q(y)} , (5.14) 

(iii) ± Hmt) {(q x T A l"d 2 q x ) + (8%^%)} + e^m^qj^d^)] 

x V(r) [q{y)T A l\{y)] , (5.15) 

(iv) ± (- e6 ( mf ))_jLfe*** ldpq x T A lv d a q x ] V(r) [q(y)T A ^ »g(y)] . (5.16) 

Here, sign + (— ) means that quark chirality is R (L). In the non-relativistic limit, (i) and 
(ii) vanish, since components of fi = and \x = i are cancelled with each other. For this 
calculation, we have used a NRQCD result, d°q ~ C((m^) 3//2 ) (v: c-quark velocity, m c : c- 
quark mass). Actually, (iii) induces the leading effects for the potential. By taking leading 
order of v, a power counting shows 

$V™ G ( r ) = [MmgJ - e 4 (m & )) - 2(e 3 (m 3 ~J - e 3 (m & ))] 

_gj_ /-zm g \ / V{r)dl + W(r) \ ( , m 

96tt 2 V 8iV c y V V{r)& r + Vr^W ^' fc [^(r)d r fe + % k r V(r)} J' [0 - ll) 

where we use color factor (Cf/{2Nc)) from Fierz transformation. As for (iv), a, (3 must be 
space-index, so that the second derivative of space-index appears, which corresponds to D-state 
(or higher angular momentum states), so that it does not contribute the mixing between S- 
and P-states. The (iv) does not contribute the mixing between S- and P-states, too. Thus, the 
leading order of parity violating potential, which triggers the mixing between S- and P-states, 
is given by 

5V™ SY (r) = 6V^ F (r) + 5V«f(r). (5.18) 
Then, we can calculate V VciXc0 in a charmonium, and a formula of decay width is given by 



TX1X ~ 



V 2 



r(x c o^vr7r). (5.19) 
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A wave function of charmonium is given by ip{r) = ij)n,i,m(f) ■> where S-state ("01,0,0 (r)) an d 



P-state (ip 2 ,i,m(r)) are given by ipi, fi(r) = 4^ (i) 3/2 e «, ^2,i,i(r) = g^F Q)^ ^ e 2 " sinfle** 



• n 3/2 



^2,i,o(r) = (i) 3/2 2 r a CO s6», $2,1,-1 (f) = (J)^ a« sin6>e respectively Here a 



1x3/2 



8v^ 



is a Bohr-radius of the meson, and a Vc (a Xc0 ) stands for that of r\ c (xco)- Therefore, V^^J is 
finally obtined as 



ySUSY _ y(4F) , y( M G) 

»7c,Xc0 »7c,Xc0 ~ f?c,XcO' 



VcXcO 



192tt 2 8m c iV c 



(c,5) 
LL 



c 



(c,c)n 



C F (L> 



(6,5) 
LL 



(c,c)n 



4\/2 V a r), 



3/2 



(5.20) 

5/2 



"XcO 



(5.2i: 



V t q So = [M"**) - e 3KJ) - 2(e 2 KJ - e 2 KJ)] 



x 



967T 2 



-im r Cf 



IN, 



c 



( 1 



4v/2 (' 



3/2 



5/2 



2a Xc0 a ric 

2XcO + On, 



Note that m = ±1 of P-state do not contribute to the mixing, since 
zero. 



(5.22) 

integration becomes 



6 Bounds for left-right non-degeneracy of squark masses 

We are in a stage to investigate bounds for left-right non- degeneracy of squark masses. At 
first, we investigate bounds for c by use of calculation tools in this paper. Next, we estimate 
bounds for u and d by use of a similar technique in Ref.|14j. And finally, we comment on 
bounds for other sfermions. 



6.1 Bound for c 

Let us investigate the left-right non- degeneracy bound for the masses of cl and cr by use of 
the calculation method shown above. For a charmonium, we focus on r] c , whose decay has 
upper bounds of P and CP violations as [15] 

Br{j] c -»■ 7r + 7r~) < 6.0 x 1(T 4 , Br{j] c ->■ vrV) < 4.0 x 10~ 4 . (6.1) 

Note again that r] c can not decay to nn until it picks up parity violation. On the other hand, 
a branching ratio of Xco tttt is 

Br( Xc0 vttt) = (8.4 ± 0.4) x 10~ 3 . (6.2) 
A branching ratio of r] c — > nir in a direct parity violation from Eq. ( 15.191) is given by 

\F s (m )\ 2 \ib (0)\ 2 

Br( Vc -)■ 7T7r) dir . =| A uc + A dc + B cu + B cd | 2 \^ r , (6-3) 



13 



where T Vc is the total decay width of r\ c . Here we take a scalar form factor of pion F s by an 
input parameter as F s (m 2 c ) = 1, 0.1, 0.001, since its theoretical estimation is difficult above 
1 GeV. On the other hand, the indirect parity violation in rj c — >■ tttt suggests 



Br(r) c 



tttt hndir. 



VcXcO 1 VcXcO 



m 



Vc 



m 



X.cl) 



Br(x 



7T7T 



(6.4) 



where is the SM background induced from a Z-boson exchange. It gives an additional 

effect V(r) = (a/r) exp(— mzv)) in Eq. (l5.11|) . which is shown as 



SV. 



EW/ 







COS 2 



2 
3 



- - sin 2 9 



w 



lm r N 



AV(r)dl 
4^;V(r) 4ze ij ' fc §~ r fc V(r) 



(6.5) 



in a basis of (S-state, P-state) with N = 2 and = 3/2. Then, by using ip^ and V'xco' ^ ne 



SM contribution of Tj c -Xco mixing is given by 

2 



V. 



EW 

VcXcO 



9 



cos 2 6 W 



sin 2 6 



iC F 



8m c N AV2 



a XcO a Vc 



3/2 



^ a Xc0 a Vc 



a 2 X cO 



2a XcQ a r)c m z 



(6.6) 



and the branching ratio is given by Br(r] c — > titt] 



SM 



where we take E, 



yBW 



'XcO 



B r (x 



c0 



7T7T 



3.9x10 



-17 



and E Xc0 



m 



XcO- 



In Figs. 3 and 4, the branching ratios of r\ — > nir from direct and indirect parity violation 
effects are plotted, respectively, where horizontal axis is a magnitude of (m 2 L — m 2 R )/m^. 
Note that the branching ratio from indirect parity violation is larger than that from direct 
parity violation. Unfortunately, we can show that the SUSY parity violating effect is smaller 
than the experimental bound of Eq. (16. ip in the parameter region, and it is difficult to obtain 
the non-degeneracy bound between rric L and vrtc R - Figures 5 and 6 show a case that g and 
£r are degenerate around 400 GeV in mass. The magnitude of the horizontal axis is varied 



from (mk r 



m% 



J/m], = 7.8, which is taken to be consistent with LHC data. Notice that the 



branching ratio becomes larger than that in Figs. 3 and 4, however, the experimental bound 
is also much higher, and we can not obtain the bounds. 



6.2 Bounds for u and d 

The left-right non-degeneracy bounds for u and d was studied by use of nuclear parity violation 
in Ref.[H]. Where they compared coefficients of (quark level) meson- nucleon couplings in the 
SM with those in the SUSY. However, studied parameter region was m? mj < O^Gp 1 ), 
which is already experimentally excluded, so that we investigate the left-right non- degeneracy 
bound in a wider parameter region, besides, without approximations used in Ref.[H]. 
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\F 3 = 0.01 



4.X10 -17 



10 



SM 



rn~ — m% 



Figure 3: Branching ratios of rj — > 7T7t from di- 
rect parity violation with = 700 GeV, rriu R = 
1000 GeV, mu L = 2000 GeV, m rf - R = HOOGeV, 
m ,- = 2100 GeV, and m £R = 1200 GeV. 



Figure 4: Branching ratios of r\ — > tttt from indirect 
parity violation with m§ = 700 GeV, and mg H = 
1200 GeV. 



Br(r; c -> 7r7r) dil 



F s = 1 



SM 



\F" =0.1 



F s = 0.01 




Figure 5: Branching ratios of r\ — s- tttt from indi- 
rect parity violation with wig = 400 GeV, mj s = 



410 GeV, mar 



1000 GeV, m,- 



420GcV, 



mg L = 1100 GeV, and m~ CR = 430 GeV. 



Figure 6: Branching ratios of r\ — > tttt from indirect 
parity violation with m g = 400 GeV, and mg fl = 
410 GeV. 
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We use 7r, ui, p and nucleon couplings for the meson-nucleon coupling. The notation of our 
dimension six operators corresponds to 

G 2 (mq,mg>) , , Gi(m g -,mg/) 

5—2 = fi{mq,mg>), — -2 = /2 m 9 -,m ? , (6.7) 

3m~ 3m~ 

in Ref . |14] . where we neglect flavor mixings and squark left-right mixings (A-terms). On the 
other hand, coefficient of q-q-G vertex is written by 



C(m 2 - /m~) 43m? — 144m~m? + 153m~m| — 6 (2m? — 9m?m| + 6m?) log f 3 ) ~~ ^2m 



2 

m^ 



l g 54 (m~ — m|) 

(6.8) 

in a massless approximation of u- and <i-quarks. By using above equations, we calculate 
bounds from the SM as 

(i) |C p (vr) + C!(7r)|<|C SEW (vr)|, 

(ii) |(7» + C»|<|C SEW (a;)|, 
R |^(p) + ^(p)|<|C SEW (p)|, 

which are shown in Figs. [7J El and [91 C(7r), C(lo), and C(p) are parity violating effects 
(coupling) from 7T-, a;-, and p- nucleon interactions, respectively. Indcies p and 6 stand for 
penguin and box diagram contributions, respectively. Index SEW means an effect from the 
SM electroweak interactions [Hj as |C sew (tt) | = 8.5 x 10~ 7 , |C SEW (a;)| = 4.5 x 10" 6 , and 
|C SEW (p)| = 6.2 x 1CT 7 . The factor c(m^) is defined by c{rriq) = C(m^/m^)/m^, and then 

A a 2 

C P 0) = H m u R ) - c{ma L ) - c(m^) + c(mjj] , (6.9) 

C,P(W) = \ 24 P ^ (mSfl) ~ C(msi) + c(m <^ ~ C(m ^^ ' (6 ' 10) 
C p (p) = ~p [c(m Sfl ) - c(m fi J + c(mjj - c(mjj] , (6.11) 

00 = -^P [/i0"fii,™jj - h( m u R ^ m d L ) - fzijnu^m^) + / 2 (m GjJ , mjj] , (6.12) 
3a 2 /2 8 \ 

C+M = -- ^ ( - + — J p [2/1(771^,771^) - 2/i(m fiil ,m d -J - / 2 (m Si ,mj L ) + 2/ 2 (m fiK ,m^) 

- fi( m d L ' m d L ) - h{ m d^ m d R ) ~ fi( m u L ' m u L ) ~ /i(m fiH ,mc H ) 

+/a ("ij £ , ™d £ ) + /a 0"^ , ) + /a K ,m« t ) + /a K 8 , ™ 5h )] , (6.13) 
C+(P) = [ 2 /i( m « i; m J L ) - 2/1(^,771^) - / 2 (m Si ,m d - L ) + 2/ 2 (m fiii ,m d -J 

- fi( m d L , m d L ) ~ hi. m d^ m d R ) ~ h{.mu L ,rnu L ) - /i(m^,m^) 

+/2(m d ~ L ,m d ~J + f2{m SiR ,m SiR ) + / 2 (m« £ ,m fi J + / 2 (m fij{ , m fi J] , (6.14) 
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where we take p ~ vlO. 

In Figs. dEl and[9l we take sample points which are not excluded by ATLAS experiment [3]. 
Under m g = 700 GeV, = 1000 GeV, m Si = 2000 GeV, and mj = 1100 GeV, we change a 
value of (m| — ra| )/^| from 6.6 for the consistent with the experimental data. Unfortunately, 
in this parameter space, u and d are too heavy to obtain bounds for degeneracies between 
rriu L and m^, or, and raj ■ On the other hand, when gluino and squarks degenerate 
within 30 GeV, 7T-, w-, and p-nucleon couplings are shown inflO l [TT ] and [121 respectively. The 
magnitude of (m 2 ~ — m 2 ~ )/m~ is varied from 5.2 for the consistency with the LHC data. In 

this parameter space, u and d are again too heavy to obtain the bounds. The branching ratio 
is small because SUSY effects always have a loop factor, and it is the reason why there are 
the asymptotic values in Figs. 7~12. 



SM 



10-9 - 



SUSY 



u-N\ 



SM 



SUSY 



12 

mi — m\ 



12 

771 j — 771 2 T 



Figure 7: tt - N coupling (\C p (ir) + C b _{^)\ < 
|C SEW 0)l) with m~ g = 700 GeV, = 1000 GeV, 
m~ dR = 1100 GeV, and m Si = 2000 GeV. The mag- 
nitude closes in 1.0 x 10~ 8 as (m 2 T — m 2 ~ — > oo. 

v d L d R " 9 



Figure 8: u - N coupling (\C p (w) + C$.(w)| < 
|C SEW H|) with m- g = 700 GeV, m 5R = 1000 GeV, 
mg R = 1100 GeV, and mz t = 2000 GeV. The mag- 
nitude closes in 4.8 x 10 -9 as (m 2 r — to 2 ; )/to1 — > oo. 



6.3 Bounds for other sfermions 

Let us comment on the bounds for left-right non-degeneracies of other sfermions. As for b, 
a total decay width of each bound state of b 6-meson has not experimentally measured yet. 
If we can know the width, the b 6-meson system can be analyzed, and a bound for a non- 
degeneracy between mi and can be calculated just as the bound between rric L and rric R 
was calculated from the charmonium. We will calculate the bounds by just replacing t] c — > rjb 
(r] b : 0~ + ) and Xco Xbo (x&o : ++ ). We hope our method is useful to give a bound between 
m~ b and m~ b from a future experiments of 5-physics. 
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Figure 9: p - N coupling (\CP( P ) + C\{p)\ < 
|C SEW (p)|) with m g = 700 GcV, m iR = 1000 GcV, 
m~ dR = 1100 GcV, and m 5i = 2000 GeV. The mag- 
nitude closes in 1.1 x 10~ 9 as (m% — m 2 ~ )/to~ — > oo. 

v d L d R " 9 



Figure 10: tt - N coupling (|C p (7r) + C b _(ir)\ < 
|C SEW (tt)|) with Trig = 400 GeV, m fiR = 410 GcV, 
mj R = 420 GcV, and m Si = 1000 GcV. The magni- 
tude closes in 3.8 x 10 -8 as (m 2 r — m 2 7 )/m? — > oo. 
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Figure 11: a; - JV coupling (\C p (u) 
|C SEW (w)|) with m g = 400 GeV, m &R 



C b + (co)\ < 
- 410 GcV, 



420 GeV, and m iL = 1000 GeV. The magni- 



tude closes in 1.3 x 10 8 as ( 



2 



Figure 12: p - N coupling (\C p {p) + C\{p)\ < 
|C SEW (p)|) with nig = 400 GeV, m iR = 410 GeV, 



— 420 GcV, and ma L 
tude closes in 2.6 x 10~ 8 as (r 



1000 GcV. The magni- 

- m 2 ~ )/m 2 l -> oo. 

(in" « 
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As for s, it is difficult to estimate the bound from the same method in section 6.1. The 
reason is as follows. If we include a mixing between d- and s-quarks through the Cabbibo 
angle, this effect is too small to induce the bound between m~ SL and m~ SR because Figs. [7J 
and [9] can not give bounds for u, d, too. On the other hand, if we take s-quark as a heavy 
quark and calculate a quarkonium in NRQCD as c-quark, we might have bounds of s for left- 
right non- degeneracy from parity violating decay mode of 77 (548). Here, 77 (548) is h which 
might have a mixing with /o(600) (0 ++ ), if parity violation exists. The decay mode of /o(600) 
is dominated by 2ir. Thus, the parity violation induces 7/(548) — > hit, whose experimental 
bounds are given by 

Br(r/ ->■ 7r + 7r~) < 1.3 x lCT 5 , Br(r/ ->■ 2tt°) < 3.5 x 10" 4 . (6.15) 

However, these state are not composed only by s-quark but also it-, <i-quarks, so that the valid 
estimation is difficult. Also we should remind that mass of s-quark is about ten times smaller 
than that of c-quark which is too light to be treated in the NRQCD. 

Finally, we comment on sleptons. The lepton flavor violation (LFV) experiments require 
stringent bound of non- degeneracy among slepton flavors (generations). However, the LFV is 
suppressed when slepton masses are enough heavy even if their left- and right-handed slepton 
masses are not degenerate. That is, the left-right degeneracy is not required when sleptons 
are heavy enough. This situation is the same for squark sector as above (and also shown in 
K° — K° system, where left-right degeneracy is not required with enough heavy squarks). 

7 Summary 

The SUSY SM has a parity violation in QCD through a chiral quark-squark-gluino interactions 
with non-degenerate masses between left-handed and right-handed squarks. Since experiments 
have not show any parity violations in QCD yet, a bound for the mass degeneracy between 
left-handed and right-handed squarks should exist. We have tried to obtain this bound for 
each squark. At first, we have investigated a non- degeneracy bound between rric L and m £fl 
from experimental data of charmonium decay by use of NRQCD. Next, we have estimated 
the non-degeneracy bounds for u and d from nucleon-meson scattering data, and commented 
on other squarks. Unfortunately, our results are much below the current data, and can not 
obtain the left-right degeneracy bounds for squark masses. We hope our method is useful for 
obtaining the bounds by use of a future experimental data. 
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